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ABSTRACT. An exact equation is obtained that relates the products of two-point differences of fluid 



. velocity and those differences with the difference of pressure gradient and other quantities. The averages 

' of such products are structure functions. Equations that follow from the Navier-Stokes equation and 

incomprcssibility but with no other approximations are called "exact" here. Exact equations for structure 
functions are obtained, as is an exact incompressibility condition on the second-order velocity structure 
function. Ensemble, temporal, and spatial averages are all considered because they produce different sta- 
tistical equations and because they respectively apply to theoretical purposes, experiment, and numerical 
simulation of turbulence; those applications are addressed herein. The midpoint and the difference of the 
two points at which the hydrodynamic quantities are obtained are X and r; t is time. The equations 
are organized in a revealing way by use of X, r, t as independent variables. Dependences on X and on 
the orientation of r and on t fade as the asymptotic statistical states of local homogeneity, local isotropy, 
and local stationarity, respectively, are approached. The exact equations are thus applicable to study of 
the approach toward those asymptotic states. Exact equations obtained by averaging over a sphere in 
r-space have a particularly simple form. The case of a simulation that has periodic boundary conditions 
leads to particularly simple equations. A new definition of local homogeneity is contrasted with previous 
definitions. The approach toward the asymptotic state of local homogeneity is studied by using scale anal- 
ysis to determine the required approximations and the approximate equations pertaining to experiments 
and simulations of the small-scale structure of high-Reynolds-number turbulence, but without invoking 
local isotropy. Those equations differ from equations for homogeneous turbulence. The traces of both 
exact and approximate equations have particularly simple forms; in particular, the energy dissipation 
rate appears in the exact trace equation even without averaging, whereas in previous formulations the 
energy dissipation rate appears after averaging and use of local isotropy. The trace mitigates the effect 
of anisotropy in the equations, thereby revealing that the trace of the third-order structure function is 
expected to be superior for quantifying asymptotic scaling laws. 
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The dynamic theory of the local structure of turbulence is so named by Monin and Yaglom (1975) (their Sec. 
22) to mean the derivation and investigation of equations for structure functions by use of the Navier-Stokes equation. 

O \ The structure functions are averages of differences of basic hydrodynamic quantities such as velocity and pressure 
gradient. Monin and Yaglom (1975) pointed out that the dynamic theory gives important relationships between 
structure functions, and that these relationships provide important extensions of predictions based on dimensional 
i-C . analysis and flow similarity. The dynamic theory is the basis for Kolmogorov's (1941a) famous equation that relates 
second-order and third-order velocity structure functions, and is of fundamental importance in the theory of locally 

J> , homogeneous and locally isotropic turbulence (Monin and Yaglom, 1975; Batchelor, 1947; Monin, 1959; Frisch, 1995). 
The dynamic theory does not uniquely determine the structure functions; this is known as the closure problem 
(Monin and Yaglom, 1975). Experimental data have been used to evaluate the balance of Kolmogorov's equation and 
■ generalizations of it (Antonia, Chambers, and Browne, 1983; Chambers and Antonia, 1984; Lindborg, 1999; Danaila 
et ai, 1999 a,b; Antonia et ai, 2000). This report supports such experimental work, as well as more precise use of 
direct numerical simulation (DNS) by giving correct and complete equations to be used in such evaluations. 

We derive exact equations for structure functions by use of differential operator identities. By "exact" we 
mean that the equations follow from the Navier-Stokes equation and the incompressibility condition with no additional 
approximations. This meaning is emphasized because turbulence researchers consistently use "exact" when they mean 
asymptotic. Exact equations satisfy the perceived need by Yaglom (1998) for careful derivation of dynamic-theory 
equations and the perceived value placed by Sreenivasan and Antonia (1997) on aspects of turbulence that can be 
understood precisely. In Sec. 2, the equations for products of differences is developed to the greatest extent possible 
before any average is performed. This mathematical method is similar to that used in the theory of wave propagation 
in random media where the equations for wave-field products are thoroughly developed before an average is performed 
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(see Rytov et al., 1989). Our study is limited to two spatial points and a single time and to the lowest-order equation 
of the dynamic theory; that equation includes second- and third-order velocity structure functions. On the other 
hand, the mathematical method used is of wider applicability; it is not limited to just two points, a single time, the 
lowest-order equation, or to the Navier-Stokes equation alone. For instance, the method could be used to derive an 
equation involving a three-point structure function for a scalar quantity having its continuity equation coupled to 
the velocity field. Such an equation would be useful for interpreting the observed (Mydlarski and Warhaft, 1998; 
Sreenivasan, 1991) local anisotropy of scalar fields in the presence of a mean gradient of the scalar. 

The exact equations retain all of the dependence of the structure functions on r, X, and time t, where r is 
the vector spacing between two points at which the measurements are obtained within the turbulent flow and X is 
the midpoint position of these points. Previous methods (Batchelor, 1956; Lindborg, 1996; Hill, 1997) of deriving 
dynamic-theory equations neglected the dependence of statistics on X, and thereby limited the equations to the cases 
of homogeneous and locally homogeneous turbulence. To also study the approach toward local homogeneity equations 
are needed that retain X. Here, attention is given to the conditions that must be fulfilled for the X-depcndencc to 
be neglected. Previously (Hill, 1997), the approach toward local isotropy was examined, although exact equations 
were not then available. Consequently, the approach toward local isotropy is not considered here. 

A scale analysis is performed to quantify terms that are to be neglected on the basis that |r| is much less than 
a length scale that will be called the outer scale, and to deduce all other required approximations. Such scale analysis 
is presented in detail in Sec. 7.4. Our analysis determines approximations that quantify the degree to which the 
small-scale structure of turbulence depends on its large-scale structure; such analysis was called for by Yaglom (1998). 
Our analysis sets the stage for DNS and experimental studies of the approximations. 

The equations derived in Sections 2 and 3 are exact for every flow, whether laminar or turbulent. For example, 
the equations apply exactly to the edge of a jet, to a boundary layer, as well as to those experimental situations 
such as grid-generated wind-tunnel turbulence, for which local homogeneity is expected to be most accurate. The 
equations apply provided there are no forces on the fluid at the points of measurement. Forces can be applied near the 
point of measurement; for instance, the equations are exact for hot-wire anemometer supports just downstream of the 
measurement points. The equations apply for turbulence generated at places other than the points of measurement; 
examples are grid-generated turbulence measured downstream of the grid, and turbulence generated by rotating 
blades (Zocchi et al., 1994). The case of statistically homogeneous forces distributed throughout the fluid has been 
considered for the asymptotic case of isotropic turbulence by Novikov (1965) (see also Frisch, 1995). The case of 
forces at the points of measurement is considered in Appendix A. 

The ensemble average is considered first (Sec. 3.1). It has the advantage for theoretical studies that temporal 
and spatial changes can be considered because the ensemble average does not eliminate dependence on X or t. The 
temporal average is typically used with experimental data, and the spatial average is typically used for data from 
DNS. For this reason, exact equations for both temporal averaging (Sec. 3.2) and spatial averaging (Sec. 3.3) are 
also obtained. The connection between the derivations presented here and any experiment or DNS is important 
because the equations relate several statistics and therefore are most revealing when data are substituted into them. 
A recently developed experimental method (Su and Dahm, 1996) has the potential to thoroughly evaluate terms in 
the equations derived here. As shown in Sec. 3.4, the exact equations have a particularly simple form for the case of 
DNS with periodic boundary conditions. 

The equations can be evaluated with experimental or DNS data to determine the most significant terms in the 
equations for a given flow and thereby determine the effects that cause deviations from asymptotic laws. The ongoing 
interest in turbulence intermittency includes accurate evaluation of inertial-range exponents of structure functions, 
for which purpose precise definition of an observed inertial range is needed. The third-order structure function can 
serve this purpose because it has a well-known inertial-range power law and the 4/5 coefficient (Kolmogorov's (1941a) 
4/5 law) in the asymptotic limit of accurate local homogeneity and local isotropy. Deviations from the 4/5 coefficient 
are observed in experiments (Anselmet et al., 1984; Mydlarski and Warhaft, 1996, 1998; Lindborg, 1999); this casts 
doubt on the precision with which measured exponents apply to the intermittency phenomenon. The equations 
derived here, when evaluated with data, can reveal the effects contributing to the deviation from Kolmogorov's 4/5 
law. The usefulness of such evaluations is shown by Lindborg (1999); Danaila et al. (1999 a,b); and Antonia et al. 
(2000). They generalize Kolmogorov's equation by the addition of a term describing streamwise inhomogencity. To 
obtain this term from the present exact analysis, it is necessary to perform the Reynolds decomposition. The present 
analysis has the advantage that it reveals all terms that describe inhomogeneity. This is discussed in detail in Sec. 7. 
The equations derived here are obtained in the Eulerian framework, which is most useful for experimental evaluation. 

Particular attention is given to the typical experimental case that is used to investigate universality of turbulence 
statistics at small scales and large Reynolds numbers. We derive the simplification of the exact equation that applies 
approximately to such experiments. Experimental data typically have the mean velocity subtracted before structure 
functions are calculated from the velocity fluctuation. For this reason, we derive the approximate equation obeyed 
by structure functions calculated from velocity fluctuations. The Reynolds decomposition (Sec. 5) is essential for 
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this purpose. The derivation is necessarily long in Sec. 7.4, but in this case, the journey is more significant than the 
destination because all required approximations are determined en route. Local homogeneity is the most important of 
the approximations. A necessary condition for local homogeneity is given in Sec. 7.3; it is not a sufficient condition. 

The trace of the exact equation has a particularly simple form. When averaged over a sphere in r-space, and 
when the advective and time-derivative terms are neglected, this equation has the same form as Kolmogorov's (1941a) 
equation (Sec. 4.3). This is true despite the fact that the r-space sphere-averaged equation is valid even for extreme 
violations of local isotropy. 

1.1 Contrasting Definitions of Local Homogeneity 

Local homogeneity has been given various definitions by different authors. Kolmogorov (1941b) introduced a 
space-time domain that is small compared to L and T=(L/U), where L and U are "typical length and velocity for the 
flow in the whole." Kolmogorov considers the two-point differences of the velocities at spatial points in the domain; 
one point is common to all the differences. Kolmogorov (1941b) defines local homogeneity as follows: the joint 
probability distribution of the velocity differences is independent of the one common spatial point, and of the velocity 
at the one common point, and of time. Data of Praskovsky et al. (1993), Sreenivasan & Stolovitzky (1996), and 
Sreenivasan & Dhruva (1998) contradict the statistical independence of velocity difference and the velocity at either 
end point, as well as contradict the statistical independence of velocity difference and the velocity at the midpoint. 
The exception is isotropic turbulence (Sreenivasan & Dhruva, 1998) for which case local homogeneity is assured. An 
alternative possibility that is particularly relevant here is that the two-point velocity sum, u n +u' n might be statistically 
independent of velocity difference, but statements by Sreenivasan & Stolovitzky (1996) and Sreenivasan & Dhruva 
(1998) contradict that statistical independence as well; publication of supporting data would be useful. Kolmogorov's 
definition should not be used because experimental data contradict that statistical independence (Praskovsky et al, 
1993; Sreenivasan and Stolovitzky, 1996; Hill and Wilczak, 2001), as do theoretical considerations (Hill and Wilczak, 
2001). 

Monin and Yaglom (1975) define local homogeneity to mean that the joint probability distribution of the two- 
spatial-point velocity differences is unaffected by any translation of the spatial points. They do not impose a restriction 
on the translations to a spatial domain. It follows (Monin and Yaglom, 1975) that statistics composed entirely of the 
differences obey the same relationships that they do for homogeneous turbulence (namely, they arc independent of 
where they are measured), and that the mean velocity depends linearly on position. In practice, statistics of differences 
and of derivatives do depend on where they are measured except in the ideal case of homogeneous turbulence. Frisch 
(1995) gives a definition that is equivalent to that of Monin and Yaglom (1975), except that the translations are 
restricted to a domain the size of the spatial scale characteristic of the production of turbulent energy (which he calls 
the integral scale). Two-point structure- function equations of all orders contain a statistic that is the product of not 
only factors of the difference of the two velocities but also one factor of the sum of the two velocities, i.e., u n + u' n 
(Hill, 2001). Because the definitions of local homogeneity by Monin and Yaglom (1975) and Frisch (1995) involve 
only the joint probability distribution of two-point differences, it follows that those definitions are not sufficient to 
simplify structure- function equations to the same level of simplification as does homogeneity. 

The calculus of homogeneity by Batchelor (1956) is the commutation of spatial derivatives from within an 
average to outside the average where they become derivatives with respect to r, and vice versa. The calculus of local 
homogeneity by Hill (1997) is a generalization of Batchelor's calculus; specifically, local homogeneity was implemented 
by neglecting the derivative with respect to X relative to the derivative with respect to r when spatial derivatives 
were commuted with the averaging operation. That implementation is restricted to statistics that contain at least 
one difference or derivative of basic hydrodynamic quantities (such as velocity, pressure, temperature, etc.). This 
calculus differs from the aforementioned definitions of local homogeneity in that no translational invariance is required 
other than for the infinitesimal displacement in X implied by the derivative operation. In Appendix C, examples are 
given that show how this calculus produces the predictions of homogeneity for the homogeneous case. To simplify 
the structure-function equations, Hill (1997, 2001) found that it was necessary to apply that calculus to statistics of 
products containing not only at least one difference but also quantities that were not differences. 

Consider grid-generated turbulence in a wind tunnel operated with constant mean velocity. For anemometers 
fixed relative to the position of the grid, the turbulence is stationary and streamwise inhomogeneous. For simplicity, 
ignore the cross-stream inhomogeneity. For anemometers moving relative to the grid in a direction parallel to 
the streamwise direction, the turbulence is both streamwise inhomogeneous and nonstationary. It is nonstationary 
because of downstream decay of the turbulence intensity. That example raises the question as to whether or not local 
stationarity and local homogeneity should be combined into a single definition that is independent of the motion of 
the coordinate system. In this author's opinion such a combined definition is neither desirable nor practical. Thus, 
local homogeneity (or local stationarity) must be considered in a given coordinate system. 
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2. EXACT TWO-POINT EQUATIONS 

Exact equations are given here that relate two-point quantities and that are obtained from the Navier-Stokes 
equations and incompressibility. The two spatial points are denoted x and x'; they are independent variables: they 
have no relative motion; e.g., anemometers at x and x' are fixed relative to one another. To be concise, velocities are 
denoted Ui = itj(x, t), u[ = Wj(x', i), and the same notation is used for other quantities. p(x, t) is the pressure divided 
by the density (density is constant), v is kinematic viscosity, and d denotes partial differentiation with respect to its 
subscript variable. Summation is implied by repeated Roman indices; e.g., d Xrl d Xrl is the Laplacian operator. For 
brevity, define: 

dij = (ui - m-) (uj - u'j) ; (1) 
dijn = (u-i - u[) (uj - u'j) (u n - u' n ) ; (2) 
nj = (d Xi p - d x >j>') (uj - u'j) + {d Xj p - d x >.p'^j {ui - u'j) ; (3) 
e%i = (d Xn Ui) (d Xn Uj) + (dx^u'i) (d x > n u' 3 ) ; (4) 

F ijn = (m - u[) ( Uj - u'^ Un ^ Un . (5) 

We change independent variables from x and x' to the sum and difference independent variables: 

X = (x + x') /2 and r = x — x', and define r = |r| . (6) 
The derivatives and d ri are related to d Xi and d x >, by 

d Xi = d ri + ^d Xl , d x * t = -d u + ^d Xz , d x% = d Xi + d x >. , d n = i (d Xz - d x >J . (7) 

It is essential to hold fixed the correct variables for each of the above partial derivative operations. The partial 
derivative d Xi is obtained with the following variables held fixed: xj, for j =/= i, and x' and t. Likewise for d x i, , x'p for 
j 7^ i, and x and t are held fixed. For d Xi , Xj, for j ^ i, and r and t are held fixed. For d ri , rj, for j ^ i, and X 
and t are held fixed. For any functions /(x, t) and c/(x',i), (Q) gives 

d rt [/(x, t) ± g(x', t)} = d Xi [/(x, t) t 3(x', t)} /2. (8) 

For example, d ri (Uj — u'-) = d Xi (uj + u'A/2 , and d ri (Uj + u') = d Xi (uj — u'j) /2. 

Now, Tij and the trace of (g) and (Q) (i.e., r« and en) can be expressed differently. Use of (0) in (||) as well 
as in en and rearranging terms gives 

nj = -2 (p - p>) ( Sij - S y + d Xz [{p ~ p') ( Uj - u'j)) + d X j [(p - p') ( Ui - u[)] , (9) 
2ve u = 2{e + e') + 2ud Xn d x „ (p + p') , (10) 



where s,j = (d Xi Uj + d X] u t ) /2 , and e = 2vs i2 Sij ; (11) 

to obtain ( |Io| ) we used Poisson's equation d Xrl d Xn p = —d Xi Ujd Xj Ui. Incompressibility requires that the trace of s^- 
vanishes; thus, the trace of @ is 

T ii = 2d Xi [(p-p')(u i -u' i )}. (12) 
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2.1 Use of the Navier-Stokes equation 



The Navier-Stokes equation for velocity component Ui (x, t) and the incompressibility condition are 

d t Ui + d Xn (UiUn) = -d Xi p + vd Xn d Xn Ui , and d Xn u„ = 0. (13) 

By multiplying the Navier-Stokes equation for by u'j, we obtain an equation having u'jdtUi as its time-derivative 
term. We add and subtract eight such equations to obtain the equation having as its time-derivative term the 
expression UjdtUi — u'jdtUi + u'jdtu\ — Uidtu'j + u\dtUj — u^dfUj + UidtUj — Ujdtu' i = dt [(tii — u'j) (uj — u'j)] . Algebra is 
used to simplify the terms in the resultant equation, and zero is added to the equation (for convenience) in the form 
of d Xn (u'iU'jUn) + d x ' n (uiUju' n ) (which vanishes by incompressibility). We thereby obtain 



d t dij + d Xn (dijU n ) + d X ' n (diju' n ) = -ry + v {p Xn d Xn d i3 + d^d^d^) 



(14) 



Use of (0) in @, and use of the identity d Xn d Xn (fg) = fd Xn d Xn g + gd Xn d x J + 2 (d Xn f) (d Xn g) to simplify 
the terms proportional to v gives 



dtdij + dx n F ijn + d rn di 



(15) 



As a check, one sees that (|l5|) is the same as can be obtained by specializing, for the present case, equation (2.13) in 
Hill (2001). The trace of (|l5j) and substitution of @ and @ give 



dtdu + dx n Fim + dr n d lin = 2ud rn d rn du — 2(e + e') + w, 



(16) 



where w = -2d Xn [(p - p) (u n - u' n )} + ^d Xn d Xn d u - 2vd Xn d Xn (p + p') . (17) 

The limit r — * applied to ( |l6| ) recovers the definition of e i n (|ll"| ) • It is significant that e appears in the unaveraged 
exact equation (nw because e will appear in the average of (|15|) only for the locally isotropic case. 



2.2 Exact Incompressibility Relationships 



Because x and x' are independent variables, d Xi v!j = 0, and d x / Uj = 0. Then, incompressibility gives: d Xn u n = 
0, d x u' n — 0, d rn u n — 0, d Tn u' n = 0, d Xn (u n — u' n ) = 0, d Tn (u n — u' n ) = 0. The combined use of incompressibility 
and (g) gives 

Q rn [{uj - W 3 ) (u n - <)] = d Xn [(Uj + u^) (u n - <)] /2, (18) 
9 rj d rn [(uj - u'j) (u n - u' n )] = d Xj d Xn [{uj + u'j) (u n + u' n )] /4. (19) 



3. EXACT AVERAGED EQUATIONS 



3.1 Exact Equations: Ensemble Average 



The ensemble is defined as a set of similar flows. An example is a set of mechanically identical wind tunnels 
operated with the same forcing. Points x and x' are defined in each flow relative to the mechanical structures or 
relative to the corresponding locations where the flow is (or was) forced. Time t is defined for each flow from the 
start of the forcing. Thus, the space-time points (x, x',i), or equivalently (X, r, t), are in complete correspondence 
between flows in the ensemble. The ensemble average is defined at each point (X, r, t) as the arithmetical average 
over the ensemble. We denote the ensemble average by angle brackets where the subscript E is a mnemonic 

for 'ensemble.' Define the following statistics: 
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Dij (X, r, i) = (dij) E , Dij n (X, r, t) = (dij n ) E , Tij (X, r, t) = \Tij) E , 

E id (X, r, t) = (e ij ) E , W (X, r, t) = (w) E , F ljn (X, r, t) = {F ljn ) E ■ (20) 

The argument list (X, r, i) is shown above to emphasize that the average applies to the general case of nonstationary, 
inhomogeneous turbulence, and that the ensemble average does not eliminate dependence on any independent variable. 
The argument list (X, r, t) is suppressed where clarity does not suffer. Defining the symbols D^, Dij ni Ty, E^j, W, 
and Fij n causes brief notation in later sections. Because the ensemble average is a summation, it commutes with 
differential operators, and the average of (jlj) is therefore 



d t Dij + d Xn F ljn + d rn D ljn = -T tj + 2v 
The average of (0) is 



(21) 



dtDu + d Xn F Un + d rn D iin = 2ud rn d rn D u ~ 2 (e + e') E + W, (22) 



where W = ~2d Xn ((p - p) K - u' n )) E + ^d Xn d Xn D u - 2vd Xn d Xn (p + p') E . (23) 



Exact incompressibility conditions on the second-order velocity structure function are given by the average of 
and |ij) as 

d rn D jn = d Xn ( (uj + v/jj K - u' n )) E /2, (24) 
d rj d rn D jn = d Xj d Xn ((uj + u' ) (u„ + u' n )) E /4. (25) 



3.2 Exact Equations: Temporal Average 

The ensemble average used above is important because it allows us to simultaneously investigate rapid temporal 
variation that a temporal average would smooth and to investigate sharp spatial variation that a spatial average would 
smooth. It is important to consider temporal and spatial averages because they are typical of experiments and DNS, 
respectively. Of course, an ensemble average can be approximated by widely separated temporal or spatial sampling 
for stationary or homogeneous turbulence, respectively. However, nearly continuous sampling is typical. Thus, we 
represent the temporal and spatial averages by integrals, but all results are valid for the sum of discrete points as 
well. The temporal average is most meaningful when the turbulence is nearly stationary, and the spatial average is 
most meaningful for nearly homogeneous turbulence. 

Let to be the start time of the temporal average of duration T. The operator effecting the temporal average 
of any quantity Q is denoted by (o) T , which has argument list (X,r,to,T); that is, 

(Q) T = I /° + Q(X,r,t)dt. (26) 

1 Jto 

For brevity the argument list (X, r,io, T) is suppressed where clarity does not suffer. The temporal average of (|l5|) is 
(dtdij) T + d Xn (Fij n ) T + d Tn (dij rl ) T — — (ry) T 



+ 2v (d rn d rn (dij) T + 2d Xn d Xn {d i3 ) T - (ey) T ^ . (27) 



The temporal average of ( |16[ ) is 

{dtdu) T + d Xn (Fiin) T + dr n (du n ) T = 2vd rn d rn (du) T - 2 (e + e') T + (w) T , (28) 

where («;) T = -2d Xi {(p-p 1 ) (it, - u-))t + 7;9 Xn d Xn (dy) T - 2vd Xn d Xn (p + p') T ■ 
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Now, (|27|) and (|28|) are exact because they are derived from (fL3|) without approximations. They differ in form from 



fl2l| ) and (22) only in that the time derivative does not commute with the temporal average. Thus, ([27]) contains 
(dtdij) T whereas ( |f5| ) contains dtDij = dt (dij) E . 

Because the data are taken in the rest frame of the anemometers and dt is the time derivative for that reference 



frame, it follows that 

(d t d ij ) T = - f ° + d^dt = [d^ (X, r, to + T)- d^ (X, r, to)] /T. (29) 

1 J to 

This shows that it is easy to evaluate {dtdij) T using experimental data because only the first (at t = to) and last (at 
t = t + T) data in the time series are used. If [dij (X, r, t + T) — dij (X, r, t )] is bounded and its ensemble mean 
varies less rapidly than T, then we can make (dtdij) T as small as we like by allowing T to be very large. 

3.3 Exact Equations: Spatial Average 

Let the spatial average be over a region R in X-space. The spatial average of any quantity Q is denoted by 
(Q) R (r,i,R), and is defined by 

(Q)m^^J J |Q(X,r,t)dX, (30) 

where V is the volume of the space region R. For brevity, the argument list (r, t, R) is suppressed where clarity does 
not suffer. The spatial average commutes with r and t differential and integral operations and with ensemble, time, 
and r-space averages. For the divergence in X of a vector q n , the divergence theorem relates the volume average to 
the surface average; that is, 

(d Xn q n ) R = ±111 d Xn q n dX = y Q I f N n q n dS^j = y £ q n , (31) 



where S is the surface area bounding the X-space region R, dS is the differential of surface area, and N n is the unit 
vector oriented outward and normal to the surface. For brevity, the notation ^ x q n is used for the X-space surface 
average in (|3l|). 



The spatial average of (15) is 



9t (^tj)n + y y F ijn + 9r n (dij n ) R — ( T ij)i 



-i' ! dr n d rn (dij) M + ~ £ d x „dij - (e y ) R ) . ( 32) 



The spatial average of is 



dt (du) m + — j> Fun + d rn (d iin ) M = 2vd Tn d rn (du) m - 2 (e + e') R + (w) M , (33) 



where (w) m = — 

v Jx„ 

The spatial average of the exact incompressibility condition (fill) is 



-2 (p - p') (u n - u' n ) + V -d Xn dij - 2vd Xn ip + p') 



9r n (d jn ) R = 2y £. («n - u 'n) (% + u'j) , (34) 

which is, on the right-hand side, a surface flux of a quantity that depends on large-scale structures in the flow. 
Similarly, (fij|) gives 

d Tj d rn (d jn ) M = — £ [d Xj (u n + u'J] (uj + u'j) . 
Of course, (133) and (t54|) are exact. 
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3.4 Spatial Average: DNS with Periodic Boundary Conditions 



The spatial average is particularly relevant to DNS. DNS that is used to investigate turbulence at small scales 
often has periodic boundary conditions. For such DNS, consider the spatial average over the entire DNS domain. 
Contributions to ^ x q n from opposite sides of the averaging volume cancel for that case such that ^ x q n = and 
therefore (dx n qn) R = 0. In (|32|) we then have ^ x Fan = and ^ x dx n dij = 0. In ( |33"1) we have / x Fu n — and 
(u>) R = 0. In (HJ), the right-hand side vanishes. Thus, in the important DNS case described above, we have the 
significant simplification that 



d t {dij) M + d rn {d ijn ) M = - {nj) M + 2v (d r7i d rn (dij) M - (e y ) R ) , 



(35) 



d t {du) R + d rn (diin) R = 2vd rn d rn (d u ) K - 2 (e + e') R , 



(36) 



and 



Proof of ct 



follows: 



dr„ (djn) R 

Using M) in 



0. 



v [(9 Xn Ui) (dx n Uj) ~ {dx' n u[) (dx^u'j)] such that d r . (ey) H = (d r] 
gral vanishes because of the DNS periodic boundary conditions and the selected averaging volume. Thus, 



we have d r . eij = d 



(37) 

xAij: where £ tj = 
this surface inte- 



d r „ 



0. 



(38) 



No approximations have been used to obtain these equations for the DNS case considered. It seems that (p5|)- 
( |36| ) offer an ideal opportunity to evaluate the contribution of the time-derivative term dt (dij) R for freely decaying 
turbulence, as well as the contribution of the pressure term (Ty) R for anisotropic turbulence, as well as the balance 
of the off-diagonal components of (p5|). 

Because we have not introduced a force generating the turbulence and because every point in the flow enters 
into the X-space average, the DNS must be freely decaying. As shown in Appendix A, it is straightforward to include 
forces in our equations. 

Performing the r-space divergence of (^) and using (p7[)-(p8|), we have 

drjdr n (dijn) R = ~d Tj (T~ij) R ■ (39) 

This exact result is analogous to the asymptotic result in Table 3 of Hill (1997). 

We can further simplify the dissipation-rate term in (|3^). Using Taylor's series, we have e (x, t) = e (X, t) + 
r fdx n £ (X, t) + ^^f-fdx^dxpE (X, <) + •••. Clearly, a great number of terms will be needed when |x — X| is outside 
of the viscous range, but the differentiability of hydrodynamic fields guarantees convergence of the Taylor series. The 
series for e (x', t) is the same as for e (x, t) in which r/2 is replaced by — r/2, such that 



e + e' = 2e (X, t) + -r n r p d Xrl d Xp e (X, t) 



2e (X, t) + d 



-r n r p d Xp £ (X,i) 



(40) 



Only terms having even-order derivatives appear in (f4(i|). The right-most term in (|40| ) has the form dx n qn, and 
therefore vanishes when averaged in X-space over the entire DNS domain for the periodic DNS case considered. 
Substituting (ff"o| ) in (|6|) gives the term 

-2 (e + e% = -4 (s (X, t)) R = -4 (e) R (t) . (41) 

The same method applied to the right-most term in ( [§51) gives 

" ( e y> R = -^{\(dx n ui) (<9x n Uj)] x =x)R (*) 



(42) 



where the subscript x = X means that the derivatives are evaluated at the point X. Of course, none of the quantities 
in (]35|)-(|42]) depends on X because of the spatial average over X. An interesting feature of (|d]) and ( [42| ) is that their 
right-hand sides clearly do not depend on r, whereas this is not obvious in ( |35|) and (|36"|). The only dependence of 
( [4l| ) and is on t. Thus, (t) on the right-hand side of j4]])-([4^) is the entire argument list. 

Of course, these results follow from the periodic boundary conditions and the fact that the averaging volume 
is over the whole periodic structure of the DNS domain. These results follow from the symmetry of that case. 
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4. AVERAGES OVER THE r-SPACE SPHERE 
4.1 Definition of the r-Space Sphere Average and the Orientation Average 



The energy dissipation rate averaged over a sphere in r-space has been a recurrent theme in small-scale similarity 
theories since its introduction by Obukhov (1962) and Kolmogorov (1962). By averaging our equations for the trace, 
we can, for the first time, produce an exact dynamical equation containing the sphere-averaged energy dissipation 
rate. The volume average over an r-space sphere of radius rs of a quantity Q is denoted by 

(QWhcrc ^ (47r(rs) 3 /3) _1 | J Jq (X, r, t) dv. (43) 

M < r s 

The orientation average over the surface of the r-space sphere of radius rs of a vector q n (X, r, t) is denoted as follows: 

I q n =( y 4n(r s ) 2 y 1 J J ^q n (X, r, t) ds, (44) 

r— sphere 

where ds is the differential of surface area, and r n /r is the unit vector oriented outward and normal to the surface 
of the r-space sphere. Both (Q) r _ sp i iere and § r q n are functions of X, rs, and t, but the argument list (X,r<j,i) is 
suppressed. In this notation, the divergence theorem is 

(d rn q n ) r _ sphcrc = q„, (45) 

Because r, X, and t are independent variables, the r-space volume and orientation averages commute with time and 
X-space averages and with X- and i-differential operators, and, of course, with the ensemble, temporal , and spatial 

averages as well. For instance, (d t (du) M ) r _ spheTe =d t «d«) K ) r _ sphcre =((a t d ll ) r _ sphcrc ) R = d t ((^} r _ sphere ) R , etc. 



4.2 Example of an Equation Operated upon by the r-Space Sphere Average 



The r-sphere average (|4 3D can operate on the structure- function equations (|2l|), (|22j), fl27|), (|2j), (p2|), (|35. 
(|35|), d36|), (|65D, (frOD, (|66|), (|67|); indeed, it can operate on the unaveraged equations jl5| ) and (|16[) as well. These 
equations have terms of the form d Tn q n ; examples are: q n = (dij n ) m , 9 rn (du) R , Du n , (dij n ) T , d Tn {du) T , etc. By 
means of ([I||), the volume average in r-space of any term of the form d Tn q n produces the orientation average of q n 
within the subject equation. After operating on (^2|) with the volume average in r-space (fl3|), the right-most term in 

that equation contains /{e + £') r -sphcrc) > which is the same as the sphere-averaged energy dissipation rate defined 

in the third equations of both Obukhov (1962) and Kolmogorov (1962) (after multiplication by 2). 

The result of the r-space sphere average of any of our equations will be clear from operating on the simplest 
equation, namely, ([56]) for the case of periodic DNS. The average of ( |36| ) over a sphere in r-space of radius rs and 
multiplication by rsj3 and use of ( [i"l| ) gives 



3 L dt((du) r _ sphcrc ) R + j {d lvn ) K = 2vj d rn (d U ) n - -^(<£) r - sp here) K - ( 46 ) 

The terms have argument list (rs,t), but ( (£) r _ sp i lcrc ) depends only on t. Of course, none of the quantities in ( |4^ ) 

. . . I K 

depends on X because of the X-space average. Despite its simplicity, (M) has been obtained without approximations 
for the freely decaying DNS case considered; (Q) applies to inhomogeneous and anisotropic DNS having periodic 
boundary conditions. 
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4.3 Kolmogorov's Equation Derived from the Sphere- Averaged Equation 



Most readers are familiar with Kolmogorov's (1941a) famous equation that is valid for locally isotropic 
turbulence. A useful point of reference is to derive it from (f46|). An index 1 denotes projection in the di- 
rection of r and indices 2 and 3 denote orthogonal directions perpendicular to r. For locally isotropic tur- 
bulence we recall that the only nonzero components of (dij n ) R are (diu) R , (<^22i)k = (d33i) R , and of (dij) R 
are (c?n) R , and (g?22)r — (^33)r- These components depend only on r such that there is no distinction in 
an r-space sphere average between r$ and r; thus, we simplify the notation by replacing rg with r. The 
isotopic-tensor formula for (dij n ) R gives (du n ) R = {r n /r) ((dm) M + 2 (^221)^) = {r n /r) (dm) R , substitution of which 
into (Q gives § rn (du n ) M = (r„/r) (du n ) m = (r n /r)(r n /r) (d ul ) R = (dui} R . Since {d Tn r) = (r„/r), we have 
§ d rn (du) R = (r n /r) (d rn r) d r (du) R — d r (du) R . For locally stationary turbulence, which is the case considered by 



Kolmogorov (1941a), the time-derivative term in (46) is neglected; then (|46|) becomes 



4 

{dui) M = 2vd r (du) R - - (e) R r. (47) 

Alternatively, we can time avera ge (|46| ); then the time derivative can be neglected with the weaker conditions noted 
with respect to the smallness of (|2^); then 

«d«i>R>T = 2vd r «d«>R>r - \ « e >«>T r - ( 48 ) 

For simplicity of notation, continue with (f47j). To eliminate (dm)^ and ((cfeiV) from the expressions (da) R = (dn) R + 
2 (^22)r and {dm} R — (dm) R + 2 (^221)^, we use the incompressibility conditions |9 r (rfn) K + (du) s — (^22)r = 0, 



and rd r (dm} R + (dm) R — 6 (d22i) R = 0, which are valid for local isotropy (Hill, 1997). Then (|47|) becomes, after 
multiplying by 3r 3 , d r (r 4 (tiin) R ) = 6vr 3 d r [r~ 2 d r (r 3 (^ii)r)] —4 (e) R r 4 , which is then integrated from to r. After 
the term proportional to v is integrated by parts and the resultant equation is divided by r 4 we have Kolmogorov's 
equation 

4 

(dm) R = §vd r (d u ) m - - (e) R r. (49) 

Two integrations over r were required to obtain the equivalent of (^9|) in section 6 of Hill (1997), whereas one integration 
over r was required here to obtain (|49|); the reason is that the r-space sphere average replaced the first integration. 
Kolmogorov's 4/5 law, (dm) R = — ^i,£)^r, for the inertial range immediately follows from (j49|). For the viscous 

range, (dm) R can be neglected in ( f49| ) such that the known relation (e) R = (15^/2r)<9 r (^ii)r = 15^ ^(9i!«i) 2 ^ is 
obtained, where the viscous-range asymptotic formula (rfn) R = ((d Xl ui) 2 \ r 2 was used. 



5. REYNOLDS DECOMPOSITION 



The Reynolds decomposition separates any hydrodynamic variable into its mean value and fluctuation and is 
essential when considering hot-wire anemometer data. In the next section, the Reynolds decomposition is used to 
elucidate the meaning of dx n Fij n , and in Sec. 7.4 to perform the scale analysis. 

For the ensemble average, the Reynolds decomposition of m(x,t) is defined by 

Uj(x, t) = Ui(x 7 t) +Ui(x,t) , where Ui(x,t) = (ui(x,t)) E , and (ui(x,t)) E = 0, (50) 

and similarly at the point x'. For brevity, U[ = f/i(x',<), and = Ui(x',i), etc. Using (f7j), the incompressibility 
condition gives 

dx n u n — 0, d Xn U n = Q, d Xn u n — 0, d Tn u n — 0, d Tn U n — 0, d rn u n = 0, (51) 

and similarly for u' n , U' n) u' n . 

For the time average (psj), the mean velocity is Ui(x.,to,T) = (ui(x,t)) T ; as in (26) this notation emphasizes 
that the mean depends on the start, to, and duration, T, of the time average, as well as on x. The Reynolds 
decomposition is Uj(x, t) = Ui(x,to,T) + Wj(x, t, to, T), such that (ui(x.,t,to,T)) T = 0. Clearly the fluctuation, 
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Ui(x,t,to,T), also depends on to and T as well as on x and t, and ( pl| ) is valid because the time average commutes 
with spatial derivatives. 

For the space average ([50]), it follows from (B) that when the integral over X operates on a single-point quantity 
like Uj(x, t), it is an integral over x such that (|30|) produces a function only of t but not of X or r. Thus, the mean 
velocity is Ui(t,W) = (iti(x, t)) R ; as in ( |30| ) this average depends on the centroid and shape of the averaging volume, 
but this dependence is not denoted explicitly. The Reynolds decomposition is Ui(x,t) = Ui(t,M) + Mj(x, t,R), which 
gives (?/i(x, t, R)) R = 0. Clearly, ( pi] ) is valid for the space average. 

For brevity, the arguments of mean quantities are not shown in the following. 



6. MEANING OF THE TERM d x „F if 



The Reynolds decomposition ( |50| ) used in the second term of ( pl| ) (i.e., dx n Fij n ) combined with ( |5l| ) gives 



dx n Fijn — 



U' 



J n 



■J" 



(52) 



where, for brevity, we define 



(Ui - U[) , f m = ( ( 



(53) 



Note that o means that the quantity is a fluctuation, e.g. Uj, and that a statistic is calculated from fluctuations, e.g., 



T in and 



Also, D, 



%3 appears m 
Consider the first term in 



not Di 



namely | (U n + U' n ) dx n Dij- If the mean flow is spatially uniform to the 
extent that U n and U' n are equal, then \ (U n + U' n ) dx n Dij becomes the same expression that Lindborg (1999) [his 
Eq.(8)] deduced as an addition to Kolmogorov's equation. His deduction was based on Galilean invariance applied to 
a uniform mean flow. The combination of (^2|) and (|2l]) shows that both dtDij and | (U n + U' n ) dx n Dij must appear 
in the dynamical equation as was correctly deduced by Lindborg (1999) on the basis of mean-flow Galilean invariance, 
but replacing d t Dij with ^ (U n + U n ) dx n Dij, as was done by Danaila et al. (1999 a,b) on the basis of Taylor's 
hypothesis, does not preserve that invariance. Now, dx„Dij is a measure of inhomogeneity because dx n Dij is the 
rate of change of Dij (X, r, t) with respect to where the average is performed. Thus, \ (U n + U' n ) dx n Dij describes 
the effect of the fluid moving relative to the anemometers in a direction in which D^- (X, r,t) is inhomogeneous. 
Lindborg (1999) quantifies the contribution of \ (U n + U' n ) dx„Dij to Kolmogorov's (1941a) equation (Sec. 4.3) for 
several experiments and thereby shows that the contribution can be significant. 

Now, i (U n + U' n ) dx n Dij is well illustrated by the case of turbulent flow in a pipe or wind tunnel. Perform 
the X-space spatial average (|30| ) of \ (U n + U' n ) dx n Dij over a cylinder having sides parallel to the mean velocity and 
having ends perpendicular to the mean velocity. For simplicity, assume that the mean velocity is uniform over the 
ends of the cylinder so that U' n = U n = |U| s n where s n is a unit vector in the streamwise direction, which is the 
1-axis. Use of the divergence theorem (pi]) gives 



1 

V 



1 



{Un + U^Dij dX 



1 
AL 

H 

L 



iV„s„ |U| DijdS 



SnDij 'downstream 



Sfl -Dij ^upstream 



(54) 



where ^ x s n D. 



ij downstream 



and 



n-^ijupstream 



are the surface averages over just the downstream and upstream 



ends of the cylinder, respectively, and L, A, and V = AL arc the length, area of the ends, and volume of the cylinder, 
respectively. Now, (|U| /L) 1 is the mean time for the flow to pass from the upstream end of the cylinder to the 
downstream end. Thus, (p4j) is the rate of downstream decay of D,j averaged over the cylinder cross section. 

Now consider the term dx„ (^A-iTj n + AjT; n + L jn j in (|52|). From ( |53|) this term is important if there is 

strong correlation between velocity difference and velocity sum. One such case is when at least one anemometer is 
at the edge of a jet and is therefore sometimes immersed in quiescent entrained fluid and sometimes in turbulently 
agitated fluid. More generally, the second term in ( p32| ) is important for the case of large-scale structures. This 
term describes a contribution caused by inhomogeneity in the direction transverse to the mean flow direction as well 
as in the streamwise direction. Thus, this term is expected to contribute for pipe and jet flows when anemometers 
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are separated transverse to the flow. Experimental and/or numerical evaluation of these terms is needed to quantify 
their contribution to ( ^l|) for particular flows. 

On the other hand, the second term in (|52]), i.e., dx n (AiTj n + AjT in + Tij n ^j ,does not grow if \ (U n + U' n ) 

increases, as does the first term, i.e. , Un ~^ u ™ dx n Dij . Therefore, for a flow in which large-scale structures are minimized, 
such as grid-generated turbulence, and for a large enough Reynolds number such that r can be much less than the 
integral scale, the second term in (]5^) is expected to be negligible because it is dx„ operating on an average. For such 
a flow, one expects that the two-point sum, (u n + u' n ) , has a weak statistical relationship to the difference, (ui — u'^. 
The negligibility of the second term in ([52]) when ( p2[ ) is substituted in (^l|) will be considered further in Sec. 7.4. 



7. APPROXIMATE EQUATIONS PERTAINING TO EXPERIMENTS ON THE SMALL-SCALE 
STRUCTURE OF HIGH-REYNOLDS-NUMBER TURBULENCE 



We are now in a position to investigate three closely related objectives that will be considered simultaneously. 
One objective is to study the simplification of ( ^lf ) on the basis of data for the small-scale structure of high-Reynolds- 
numbcr turbulence; another is to determine the approximations required for that simplification. The third objective 
is to obtain from ( |2l| ) an equation that is closer to the measurement process of extracting a mean velocity from 
anemometry data. We use the ensemble-averaged equations because they retain both temporal and spatial variability. 
Here, we consider the approach toward local homogeneity. For this purpose, our equations that depend on the location 
of measurements, i.e., X, are needed. We also consider the approach toward local stationarity, so the dependence on t 
is needed. The restrictions required by local isotropy are not used, so dependence on the orientation of measurement, 
i.e., r/r, is retained. On the other hand, assumptions about the order of magnitude of some quantities require that 
local isotropy is not greatly violated. The data used for this investigation are given in Appendix B, which includes 
the empirically verified (Monin and Yaglom, 1975) formulas for the inertial and viscous ranges for components of Dy 
and D lm . 



7.1 Structure Functions of Fluctuations 



An experimenter usually extracts Ui from the anemometer's signal, then calculates statistics from Ui, e.g., 
= ((ui — u'i) (v,j — . Similarly define D^ n , Tij, and Eij in terms of the fluctuations of velocity and pres- 



sure. However, (u n + u' n ) /2 in (^|) cannot be replaced by (u n 
replacement would result in i<y n being defined as Vij n in (|53|). 



- u' n ) /2 without destroying the meaning of Fij n ; that 
A reasonable choice for the symbol Ffj n is 



F, 



U n + U' 



U n + U' 



D; 



Now ( pi] ) is not exactly satisfied by substitution of D^, Fi jn , 



Dij n , T„, and E tj , in place of Dy, F ijn , 



D 



T^, and Eij, nor does that substitution satisfy any equations derived from (21), Kolmogorov's equation being one 



tjn , 



such equation (see Hill, 1997). Substitution of the Reynolds decomposition of D^, Fij n , Dij n , Tij, and E^ (e.; 
D^ = A^Aj + D^, etc.) in ( pl| ) gives a complicated equation. Below, simpler approximate equations are derived by 
scale analysis and are summarized in Sec. 8. 



7.2 Experimentally Evaluatable Exact Incompressibility Conditions 



Because the approximations d Tn Di n — and d rn D; m — have an essential role in many theories, experimental 
evaluation of these approximations is desirable. However, the expressions d Tn Di n and d Tn Di n are nearly impossible 
to evaluate experimentally. Use of (||) and ( |5l| ) gives exact expressions for them that can be more readily evaluated; 
namely, 

Or n D in = ([d rn (ui - %)] (u n - u' n )) E = d Xjl {{ui + %) (u„ - u' n )) E /2, (55) 

which is similar to (Q). For the temporal average, the right- most expression in ( |55| ) requires, at most, mea- 
surements at four positions of the statistic ((ui + u' i ) (u n — u' n )) T . If, as in the case of grid-generated turbulence, 
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inhomogeneity is streamwise, then only two positions displaced in the streamwise direction suffice to determine 
®Xi {{ui + «■) (ui — u[)) E /2. The Reynolds decomposition gives d Tn D in — A n d rn Aj + d rn D in , which shows that 
evaluation of d Tn Di n only requires mean velocity measurements at several positions in addition to the previous 
evaluation of d rn Di n . 



7.3 A Necessary Condition for Local Homogeneity 



We must define several scaling parameters determined by the flow. The integral scale, as traditionally defined, 
is strictly applicable only to homogeneous turbulence; see, for example, Tennekes and Lumley (1972). Here, however, 
we are studying inhomogeneous turbulence. As an example of the difficulty of defining integral scales in general 
inhomogeneous turbulence, consider the horizontally homogeneous atmospheric surface during daytime convective 
conditions. It is difficult to imagine a useful integral scale defined using data obtained along a line from the ground 
to the upper reaches of the surface layer. However, the horizontal homogeneity and Taylor's hypothesis allow integral 
scales to be defined for all three velocity components measured at a point. Using surface-layer data, Kaimal et al. 
(1976) show that the horizontal velocity components scale with the depth of the entire boundary layer; that depth 
can be 1 to 2 km. Unlike the horizontal velocity component, the vertical velocity variance obeys Monin-Obukhov 
similarity such that its integral scale is proportional to the height above ground (Kaimal et al., 1976). For our study 
of the approach toward local homogeneity, it is necessary to define the large scale as the smallest of the integral scales 
or of the distance to boundaries. From the example of the atmospheric surface layer, that scale is the height above 
ground. Denote this chosen length scale by L and call it the outer scale. This name distinquishes it from the integral 
scale, which might not exist as traditionally defined in terms of the integral of a velocity correlation function. It is 
useful to define a velocity scale v by 

v^{{e) E Lf\ (56) 

Monin and Yaglom (1975) and Tennekes and Lumley (1972) determine that v is an estimate of the root-mean-square 
velocity, and that the mean shear is not greater than v/L. If this is not so for our chosen outer scale L, then L can be 



adjusted to make it so. From studies of nearly homogeneous turbulence, the right-hand side of (56) is proportional 
to velocity variance and the proportionality constant is independent of Reynolds number at high enough Reynolds 
numbers (Sreenivasan, 1998; Pearson, Krogstad, and van de Water, 2002). The proportionality constant is of order 
unity and depends somewhat on the large-scale structure of the flow (Sreenivasan, 1998; Pearson, Krogstad, and van 
de Water, 2002). 

We define the scale £ by 

£ = lOry; (57) 

£ is a scale typical of the energy dissipation range (Appendix B.l). Here, Kolmogorov's microscale 77, which is a scale 
typical of the viscous range, is defined by 

If the data have an inertial range, then £ is closely related to the r at which asymptotic formulas for the inertial and 
viscous ranges are equal; this is demonstrated in Appendix B.l. 

The basic tenet of local homogeneity is that as r is reduced relative to L, nonlinear randomization causes 
statistics of differences of basic hydrodynamic quantities to decrease their dependence on the large-scale Sow structure. 
For r < £ and as r is further reduced, the nonlinear randomization is increasingly opposed by the smoothing effect 
of viscosity. Therefore, £ <C L is a necessary condition for local homogeneity. For r > £, r <C L is the necessary 
condition. That is, local homogeneity applies to the asymptotic case: 

if r < £, then £ < L; if r > £, then rCL; i.e., max (r, £) < L, (58) 



where L is the outer scale. We study the approach toward local homogeneity by using fl58|) in scale analyses. We 
do so in Sec. 7.4, and find that some predictions of local homogeneity (such as dx„Fij n = and d rn Dj n = 0) do not 
follow solely on the basis of the necessary condition (^). Thus, ([58|) is not a sufficient condition for local homogeneity. 

Suppose for the moment that the turbulence under investigation is sufficiently homogenous that an integral 
scale L can be defined in terms of an integral of the velocity correlation function. The microscale Reynolds number 
(Tennekes and Lumley, 1972) R\ is well known to be related to integral scale L and r\ by L/rj cc R 3 ^ 2 (Tennekes and 
Lumley, 1972). Then, (M) gives L/i cc R 3 X /2 . Now, £ < L is a necessary condition in ( |58| ) ; so R\ ^> 1 is a necessary 
condition for local homogeneity, but it is not a sufficient condition. In a general inhomogeneous turbulence case, we 
assume that this is also true when L is the outer scale. 
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7.4 Scale Analysis 



This section uses the data given in equations (Bl) to (B13) of Appendix B.2. Those equations are distinguished 
by the prefix B. 



Now, we consider the scale analysis of ( |2l| ) 
local shear at point X by 



First, consider the Reynolds decomposition of D^. Denote the 



G ijn = d Xn Ui(K,t). 

On the basis ([58J) that max (r, I) <C L, we retain only the first two terms of the Taylor series of Ui and U[ around point X 
to obtain that Aj ~ r p Gin = rGi t i where the 1-axis is parallel to r. Therefore, d Tn Aj ~ (c? r „r p ) Gj iP = 8 pn Gi_ p = G^ n , 
which also follows from (|8|). Recall that the velocity scale v is defined such that (e) E is of order v 3 /L and a component 
of mean shear, i.e., Gj iP , is at most of order v/L. The Reynolds decomposition of Dij gives 



D, 



a, a, + d %3 ~ re;, ; + n, 



(59) 



Now Gj,i might be zero; if not, it is no greater than of order v/L. Use of (Bl) gives r G at \G a ,i/ D aa ~ (r/L) 



(£/L) 



4/3 



,4/3 



in the inertial range, and use of (BS) gives r 2 G Qj iG Q) i/ 'D aa 
(p8|), (|59"|) gives D aa ~ Z) QQ - Therefore, (Bl) and (B3) are used below for D, 



Consider the Reynolds decomposition of t 



le term d r D; 



in the viscous range. Thus, on the basis of 

as well as for D aa . 
in (|21|). Use of incompressibility (pW) gives 



d rn D ijn = A n (d^AiAj) + (d Tn Ai) D nj + (S^Aj) D m + 



A n d rn + Aid r D Ti 



Aj d Tn Di 



d r ,Di 



(60) 



The diagonal components of (g0|) can be compared with the diagonal elements of JSy, which, according to our data 



(B9), are of order (e) E . The first term in (|60|) can be approximated by r 2 G nj i (G aj iG Q;ra + G at iG a , n ), which is at 



most of order r 2 v 3 /L 3 ; this is of order (r/L) 2 relative to (s) E . Hence, when (60) is substituted in (|2l|), the first 
term in (pG) can be neglected relative to the diagonal element E aa on the basis of (58). The second, third, and 



fourth terms in (|6C|) introduce off-diagonal elements of D^ into the diagonal elements of (|2l|). Using (Bl) and the 
assumption (see Appendix B) that the off-diagonal elements of D^ are no greater than the D, 

and fourth terms in (60) can be shown to be no greater than of order (r / L) 2 ^' 



the second, third, 
relative to (e) E , and are therefore also 

neglected on the basis of (|58|). The same procedure can be used for the fifth and sixth terms in (|60|). On the other 
hand, substitution of the definition (|53] ) of T n j in j55|) and use of (|B13| ) gives d Tn D n j = dx n ^nj < v 2 /L such that the 
fifth and sixth term s in ( p0[) are much less than (r/L) (v 3 /L) and are therefore negligible compared with (e) E on the 
basis of both ( B13 ) and (p8[). Therefore, for the projection of J2l|) in an arbitrary direction a, (|5j|), incompressibility, 
and our data imply that ai<ijd rn Dij n can be replaced by aidjd rn Dij n . 



The stronger conclusion that ai<ijd rn Dij n 



Dij n can be obtained as follows. 



For an inertial range, 



the above comparison of terms with E aa is equivalent to comparison with ai<ijd rn Dij n because of (B5) and (|Bjj). 



Therefore, the above scale analysis is sufficient to state that for the ine rtial range aidjd rn Dij n ~ diajd r „ Du n . We 



need only extend this result to the viscous range as follows. Use of (B7) shows that the first term in (ppj) is of order 



(£/L) relative to cndjd rn Dij n in the viscous range, and that the second through sixth terms in ( p0|) are of order 



(£/L) 2 ^ 3 relative to didjd r 



Dij n . Therefore, (|58|), incompressibility, and the empirical formulas (|B5|) -(p9[) give 



Q>iQ>jO rn Dij n — QidjO rn Dij n . 



(61) 



The significance of there being a projection in an arbitrary direction a within (|6l|), is that empirical evidence is lacking 
for the off-diagonal components of Dij n . 

We are now ready to consider in more detail the second term in (52), namely dx n ( AiTj n + A^Tm + Tij n ^j . 

It is assumed that our data are chosen to mitigate large-scale structures such that ( |B13| ) is true. One part of the 

second term in ( |52"| ) is dx n AiTj n = (dx n Ai)Tj n + Aj (dx n Fj n ^J. Now, (dx n Ai) ~ rdx n GiX, this is at most 

of order (r/L) (v/L). Therefore, the ratio (dx n A a )T an / E aa is at most of order (r/L) (T an /v 2 ^J, which is very 

small compared to unity on the basis of (|58| ) and ( B13 ). Similarly, A a (dx n ^mij /E aa is of order (r/L) (r a n/v 2 ^J. 

Another part of the second term in (|52|) is dx„^ij n - The ratio ( 9v„ r acm ) / E aa is at most of order T aan /v 3 , which 
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is very small because of (B13). Therefore, the entire second term in ( |52| ) is negligible compared to E aa , and therefore 
it is negligible in diagonal components of (|2l]) . Neglecting the second term in ( |52| ) and using incompressibility, in the 
diagonal components of ( pl| ) we have 



d x ,F aan (X, r, t) ~ - ([/„ + U' n ) d x ,D aa . 



(62) 



The Reynolds decomposition of fe2j) is 



dx n F aan ~ Af QQ + - (C/„ + U n ) dx n D aai 



(63) 



where M aa = i (£/„ + U' n ) r 2 d Xn (G aA G aA ) . 

There are clearly flows for which we expect that M aa is negligible; an example is freely decaying grid-generated 
turbulence in a wind tunnel for which G Qi i = 0. On the other hand, M aa might not be negligible in all cases. 
Consider that M aa is at most of order (|U| /v) (r/L) 2 relative to E aa . Although (r/L) 2 <C 1 follows from (|5"§|), 
|U| /v can be much larger than unity. Thus, M aa cannot be neglected relative to E aa on the basis of ([58]); the same 
is true for i (U n + U' n ) dx n D aa because it is also proportional to |U| /v. We assume that the mean flow does not 



have an abrupt change near the positions of the anemometers. Then, use of (Bl) and (B3) shows that r 2 G ctj iG Q , j i 



is of order (r/L) 4 ^ 3 and (l/L) 4 ^ 3 relative to D aa in the incrtial and viscous ranges, respectively. However, it is not 
clear on this basis that we can neglect M aa relative to \ (U n + U' n ) &x n D aa because what is needed in ( |63| ) is the 

streamwise rate of change, i.e., (U n + U' n ) dx n , operating on both r 2 G a ,\G a ,\ and D aa . Consequently, we will not 
further simplify (|62]). 

Now consider the term dtDij in (pi]). Recall that the positions of the anemometers, namely x and x', are 
held fixed for the time-derivative operation dt- Thus, the meaning of dtDij is the time rate of change of Dij in the 
anemometer's rest frame. The sum of d t D aa and \ (U n + U' n ) dx„D aa [see (|6^)] is the time rate of change of D aa 
in the reference frame moving with velocity (U + U') /2; that is, moving with the fluid in the sense of moving with 
the local and momentary ensemble-averaged velocity. Now ( |2~L| ) is exact and therefore describes cases that include 
rapid changes of mean conditions in the rest frame of the anemometers. However, assume that the experimenter 
has chosen a case for which mean conditions are nearly constant in the anemometer's rest frame; examples include 
fixed anemometer positions in a wind tunnel, pipe, or jet for constant mean flow, or freely decaying DNS. From the 
Reynolds decomposition ([39]) we have 

d t D aa = r 2 d t {G a ,iG a ,x) + d t D aa . (64) 

For example, consider the case of turbulence that is freely decaying in the anemometer's rest frame, or freely decaying 
DNS. In this case, r 2 d t (G a .iG a ^) is at most of order (r/L) relative to E aa , whereas for the inertial and viscous 

•"•"v 2/3 2/3 2 

ranges dtD aa is at most of orders (r/L) 1 and (i/L) 1 (r/L) relative to E aa . For this case, d t D aa can be neglected 
in (pip. More generally, d t D aa is negligible because the experimenter chooses not to move the anemometers rapidly 
through positions where mean conditions differ greatly. Given the opposite choice, dtD aa would not be negligible; it 
would be of order (r/L) 2 (|V| /v) relative to E aa , where |V| is the speed of the anemometers relative to the large-scale 
inhomogeneous structures of the mean flow. Although (r/L) is small compared with unity, (|V| /v) can be made 
large by increasing the speed of the anemometers relative to the mean-flow structure. Thus, the term dtD aa cannot 
be neglected from ( |2l| ) solely on the basis of ( p8[ ) for the same reason that applies to dx n Fij n . We do neglect dtD aa 
on the basis of the choice mentioned above. 

Reconsider the term M aa in (|63|) together with r 2 d t (G a ,iG a ^), which appears in (64). Their sum, i.e., 



[dt + h (U n + U' n ) dx n ] r 2 (G a .iG a ^), is the temporal rate of change following the mean flow of r 2 (G a ^G a .i). This 
might not be negligible for some flows, such as a contraction in a wind tunnel or an expanding round jet, even though 
dt (G a .iG a .i) might be zero. This helps illustrate that M aa might not be negligible. 

Now consider the term proportional to v in (|2l|). The term jdx n dx n L>ij is of order (r/L) 2 relative to 
dr n dr n Dij, and is negligible. The Laplacian operating on ( [59]) gives d Tn dr n Dij ~ QGi t \Gj^\ -I- d rn d rn Dij. Now, 
2v (6Gi i iG J -.i) is at most of order 2v(v/L) 2 , which is of orders (r/L) 4 ^ 3 and (l/L) 4 ^ 3 relative to (B6) and 



respectively. Therefore, (|5|) and (B6) and (BJ) give ai&j (d rn d rn Dij + \dx n dx n Dij) ~ aidjdr^drjD^. 

In the Reynolds decomposition of E aa the terms that depend on mean velocity are of the order of an inverse 
Reynolds number (vL/v) 1 « 1 relative to (e) E . Thus, ( [f39| ) gives E aa ~ E aa . By the same method, use of the 



definition ( p"l"| ) of e gives (e) E ~ (z)e- That is, the mean velocity produces negligible viscous dissipation. 
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In the average of (]l0|), consider the term 2udx„dx„ (p + p')e> which also appears in (|23|). Excluding the case 
of nearby bodies in the flow that can cause sharp spatial variation of pressure, the mean pressure gradient scales 
with v and L. Then, the term 2vdx n dx n (p + p')e ^ s °f order (vL/v) 1 relative to (e) E , and is thus negligible. 



The Taylor series expansion (40) shows that {e + e') E ~ 2(e) E , where the neglected terms are at most of order 



{r/L) 2 relative to (s) E and are therefore negligible on the basis of (p3g|). Then, the average of ( |l0| ) gives the trace: 



) E = E ii ~E ii ~A(e) E ~±(e) E 



Finally, consider the Reynolds decomposition of Ty. Denote the mean pressure gradient at point X by 
II n = dx n (p{X,t)) E . The Reynolds decomposition of the term — 2{{p — p') (sy — s 'ij)) E hi ^ij [ see @] gives 
a mean-gradients term that is approximated by —r q r n Tl n dxq (Gj,i + Gij) = — r 2 Hidx 1 (Gj,i + Recall that 

n„ scales with v and L. Then, — r 2 H idx, {Gj,i + Gij) 1S of order (r/L) relative to E aa , such that this term is 
negligible in (pl|). Using ( |B11 ) and (B12) for the diagonal components of — 2((p — p 1 ) (s^ — s!y)) this term is 



seen to be negligible compared to E aa for r within the inertial range through the viscous range. The Reynolds 
decomposition gives dxt ({p — p') (uj — u 'j)) E — r qf n dxi (N- n Gj,q) + dx t {(p — P 1 ) (uj — u'j)) E . Since II„ scales with 
v and L, the term r q r n dxi (n n Gj iq ) = r 2 dxi (QiG j,i) is at most of order (r/L) 2 relative to E aa , and this term is 
therefore negligible in (pl|). On the basis of ( B1C ) and the neglect of —2 ((p — p 1 ) (~s a a — s^ Q; ))_E' we a l so neglect 



dx a ((p — p 1 ) (u a — u' a )) E . Taken together, these approximations show that T aa is negligible in (21). On the other- 
hand, mean pressure gradient can be large in the presence of bodies in the flow; a contraction of a wind tunnel is an 
example. Thus, like dx n Fij n , terms containing the mean pressure gradient cannot be neglected on the basis of ( |58| ) 
alone. In effect, we have assumed that there are no bodies strongly affecting the local turbulent flow. For this case, 
T aa is negligible in the diagonal elements of (^l|). 

The results of the above scale analysis are summarized in the following three sections. 

8. APPROXIMATE EQUATIONS 
8.1 Ensemble Average: Approximate Equations 

Given the experimental case discussed above and quantified in Appendix B, the diagonal elements of ( pl| ) 
projected in arbitrary directions a give the approximate equation 



- (U n + U' n ) d Xn Dij + d rn D ijn = 2vd rn d rn Dij - 



(65) 



As examples, the direction a can be chosen to be in the direction of some large-scale flow symmetry, such as streamwise 
or cross stream, etc., or in a direction defined by the separation of anemometers, such as r or perpendicular to r. 
The appearance of , rather than Dij , in the left-most term in (^5|) indicates that both terms in (^) are included. 
The trace of (^Tj) becomes 



(U n + U' n ) d Xn D u + d rn D iin = 2ud rn d rn D u - 4 (e) E . (66) 



As shown above, derivation of (|65|) and (|66j) from the exact equation ( J2l| ) requires more than just (p8[). A further 
requirement is that the experimenter avoids cases having large spatial and temporal variation of the mean flow. 
Of course, that choice improves the accuracy of local homogeneity for fixed values of [max (r, I) /L\. Additional 
requirements are approximations (B12) and ( B10| ), and that the inverse Reynolds number (vL/v) is very small. 



In general, those conditions are typical of an experimental situation that is sought for the study of the universality 
of turbulence statistics at small scales. Most experiments use Taylor's hypothesis to estimate spatial statistics from 
temporal statistics, for which purpose |U| /v must be large. For this reason, the left-most term is not neglected in 
(|65|), nor in (|66|). 

Of course, ([35|) contains no information about the off-diagonal elements of (^l[). We cannot evaluate those 
off-diagonal elements because we lack the necessary data. Clearly, DNS or a very complete experiment (e.g., as in Su 
and Dahm, 1996)) could be used to quantify those off-diagonal elements. The off-diagonal elements of (^Tj) describe 
quantities that approach zero as local isotropy becomes accurate. 
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8.2 Temporal Average: Approximate Equations 



Using ([29J), we noted the case for which (dtdij) T can be made as small as desired by use of a long averaging 
duration. This case is typical of experimental work for which the temporal average is also typical. Assume that this 
is the case such that in ( |27| ) ajdj (dtdjj) T can be neglected, and, in the case of ( |28| ) that (dtdu) T can be neglected. 
On the other hand, recall from ( p9[ ) that it is easy to evaluate a^q, (dtdg)^ from by use of experimental data. The 
Reynolds decomposition and the approximations that lead from (|2l]) and to ( |65| ) and (|6|) also apply to ( |27| ) and 
(P8|); we immediately obtain 



{u n + u' n )d Xn (d t] ) T + d rn (d 



2vd r d r (d 



it 



(67) 



K)9x n (d.. 



Or 



2vd r d r 



(68) 



where, as before, the caret over the averaged quantity means that the quantity is calculated from fluctuations. 
These equations relate the statistics that experimenters (e.g., Antonia, Chambers, and Browne, 1983; Chambers 
and Antonia, 1984; Danaila et ai, 1999 a,b) calculate from data. As shown in Sec. 5, the mean quantities, i.e., 
U n (x,tr,,T) = (it„(x, i)) T , in the definition of the Reynolds decomposition ( |50"|) are now time averages rather than 



t)) T = 0, etc. Except for replacing the ensemble average with the time average, 
and (f56|). However, the statistics in ( |65|) and (66) can have dependence on t, 
V depend on only the time of the start of the temporal average (i.e., to) and the 
duration of the average (T); in addition to which the dependence on start time and duration must be slight because 
of the neglect of (d t dij) T . 



ensemble averages such that (z 
( |67| ) and (^38j) are the same as 
whereas the statistics in (EjtT) and (I 



8.3 Spatial Average: Approximate Equations 



Now consider spatial averaging. Given the approximations that lead from (|T]) and (|2^) to ( |65| ) and (^6|) , $32} ) 
and (E3) become 



8t(d 



r + 2V J x 



(U n + U'n) dij + d rn (dijn) = 2vd Tn 8r n ( dij ) - (%) 



(69) 



dtid 



S_ 

2V 



{U n + U' n )d l% +d rn (d u 



= 2vd r d r 



-4^ 



(70) 



As shown in Sec. 5, the mean quantities, U n (t) = (m„(x, i)) R , in the definition of the Reynolds decomposition ( |50| ) are 
now space averages rather than ensemble averages such that (ui(x,t)) ffi = 0, etc. As in the previous case, the caret 

above a quantity designates that it is calculated from velocity fluctuations. The time-derivative terms dt (^dij^ and 

dt (dii) have been retained in (^9|) and (|7(]) because they are more significant than the advective term for the case 

of freely decaying DNS. Another example is the forced DNS flow of Borue and Orszag (1996), because it exhibits 
temporal variation of total mean-squared vorticity by a factor of 2. It seems prudent to retain the time derivatives. 
For DNS data, the advective term in both ([39]) and (|7(j) is seldom important. Consider the DNS flow of Borue and 
Orszag (1996), for which |U| jv was at most about 2. Then, on the basis of the scale analysis [see below (p3)], the 
advective term is negligible on the basis of (|5^). In (^) there is no information on the off-diagonal components 
because the approximations apply only to the diagonal components. 
Also, (p3) and (|7(]) become 



d t ( d, 



Or 



Hjn 



= 2vd r d r 



dtid 



= 2vd r d r 



(71) 
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9. DISCUSSION 



Given data for which local homogeneity and/or local isotropy are approximate, it seems that ( |22| ) is closer to 
that asymptotic case than is (pll), and therefore, that data for the trace Du n will more accurately show the asymptotic 
inertial-range power law than does Dm. The reason is as follows. For the approach toward local isotropy in 
homogeneous turbulence, the anisotropy quantified by nonzero values of is balanced by that from the term d rn Dij n 
in ( pl| ) (Hill, 1997). The trace of vanishes exactly for the homogeneous case because dx t {{p — p') {ui — «■)) E = 
for homogeneous turbulence and because —2 ((p — p') (su — s' u )) E = on the basis of incomprcssibility (su = 0). 
Then, d rrl D iin must balance less anisotropy in ( |22| ) than does d rrl Dij n in (pl[). For inhomogeneous turbulence, the 
nonvanishing part of the trace, namely Tu — 2dxi ((p ~ p') {ui — u' i )) E , is expected to approach zero rapidly as 



r 



decreases for two reasons. First, ((p — p') (ui — u^)) E vanishes on the basis of local isotropy. Second, the operator 
dxi causes dx t ((p — p') (uj — u 'i)) e ^° van i sn 011 the basis of local homogeneity. The right-most two terms in ( |2^ ) 
contain the operator dx„ dx n , which causes these terms in W to vanish rapidly on the basis of local homogeneity. 
Thus, all terms in W are negligible for locally homogeneous turbulence. By performing the trace it appears that 
anisotropy has been significantly reduced in ( |22] ) relative to in (|2l|). It follows that the trace, d rTl Dn n , is affected less 
by anisotropy than is d Tn Dij n , and therefore, that Du n is less affected by anisotropy than is Dij n . This hypothesis 
should be checked by comparison with DNS. Evaluation of all terms in (^||) and (|2l]) are the basis for such an 
investigation. We therefore expect that inertial-range power-law scaling would be more evident in Dm than in Dm- 
Of course, performing the trace requires that all three components of velocity be measured at both x and x'. 

To determine scaling properties of the third-order structure function, past theory has used the isotropic-tensor 
formula to produce a differential equation having the operator d r and integration of that equation (as done in Sec. 4.3). 
However, one can use an equation like (|2^) without an assumption about the symmetry properties (e.g., isotropic) of 
the structure functions by means of the sphere average in r-space, as implemented in Sec. 4.2. Evaluating resultant 
terms in the r-space sphere-averaged equation implies a tedious experimental procedure if wire anemometers are used. 
On the other hand, both DNS and the experimental method of Su and Dahm (1996) are suited to such evaluation. 
In effect, the r-space sphere average solves the equation by producing the orientation-averaged third-order structure 
function. It would seem that the orientation average mitigates anisotropy effects. Thus, the orientation average 
of the trace of the third-order structure function, namely § r Dm, is expected to best exhibit properties of locally 
isotropic turbulence, such as the inertial-range power law with the 4/3 coefficient that appears in (fl6"l). 

Lindborg (1999) estimates the contribution of i (U n + U' n ) dx n D aa (for the case U n = U' n ) to experimental 

measurements of (dm^j for grid, jet, and wake turbulence of moderate Reynolds number, and Danaila et al. 

(1999 a,b) do so for grid turbulence at R\ = 66, 99, and 448. They show that the term i (U n + U' n ) dx„D aa 
accounts for much of the observed deviation of the data from Kolmogorov's equation; Kolmogorov's equation is 

= &ud r (dii^J — | (s)t r - I n the casc of Danaila et al. (1999a), one must keep in mind that their estimation 
method reduces their equation to 2 (wi / = hm D\\ in the energy-containing range such that the balance of the 

\ / 1 — > oo 

equation is not tested in the energy-containing range. 



10. CONCLUSION 



The mathematical method of deriving exact structure-function equations from the Navier-Stokes equation is 
developed. The basic tools are the change of variables (j^) and the derivative identities (Q) and (||) and algebra. 
Manipulations are performed to the greatest extent possible (in Sec. 2) before an average is performed. Then, 
exactly defined ensemble, time, and spatial averages are used. DNS makes study of exact structure-function equations 
practical. Also, experimental methods exist (Su and Dahm, 1995) that can completel y e valuate terms in the exact 
structure-function equations. Exact incompressibility relationships, such as ( |24| ) and (pB]), are obtained. Following 
from the discussion in Sec. 9, the exact incompressibility relationship ( pi| ) will have a nonzero value at small r 
because of large-scale structures in the flow. At small r, ( |25|) is approximately the second derivative with respect to 
measurement location of the velocity covariance, and therefore clearly depends on flow inhomogeneity. 

That the exact structure-function equations are an advance can be seen from previous work. It is no longer 
necessary to derive individual terms that describe effects of inhomogeneity that arc missing from equations valid only 
for homogeneous turbulence, such as was done by Lindborg (1999). All such terms are now known. Sreenivasan 
and Dhruva (1998) note that one could determine scaling exponents with greater confidence if one has a theory 
that exhibits not only the asymptotic power law but also the trend toward the power law, and that without such 
a theory the method of computing local slopes is a "misplaced delusion." The exact equations given here are the 
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required theory for the third-order structure function, given that data must be used to evaluate the equations in a 
manner analogous to previous evaluations in Antonia, Chambers, and Browne (1983), Chambers and Antonia (1984), 
Lindborg (1999), Danaila et al. (1999 a,b), and Antonia et al. (2000). The exact dynamical equations obtained here 
are useful for studies of the approach toward local homogeneity as well as to local isotropy. Toward that end, a scale 
analysis is given in Sec. 7.4, which leads to the approximate equations in Sec. 8. The exact equations provide insight 
into the time-derivative terms, as discussed in Sec. 6. 
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Appendix A: Forced Turbulence 



The Navier-Stokes equation ( |l3| ) and the exact structure-function equations [e.g., (21)] apply to cases in which 
the turbulence is forced at places other than at the points of observation x and x', such as grid-generated turbulence, 
pipe flow, and boundary layers. Also, the Navier-Stokes equation (|l3|) and (^T|) apply to freely decaying DNS such 
as that by Boratav and Pelz (1997) and the simulation of laboratory experiments as in de Bruyn Kops and J.J. Riley 
(1998). Some DNS employ spatially distributed forces to drive the turbulence to a steady state. The Navier-Stokes 
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equation (|l3|), and the exact equations derived from it, do not apply to that case; instead, such forces must be 
introduced into ( |l3| ) and the resultant additional terms derived for the exact structure-function equations. 

If a force fi is added to the right-hand side of the Navier-Stokes equation (|ll|), then the term to be added to 
( |l5| ) is simply — ry defined in (g) with —d Xi p replaced by fi and —d x >p' by /,-. That is, the added term is 

(fi - ft) (uj - u'j) + (/, -/;) (m ~ u\) = $«,-, 

and the average of this expression must appear in our subsequent structure-function equations. Consider the case 
of the deterministic force, /j = S i2 F cos (kfXi), that was used in the DNS in Borue and Orszag (1996), where 
we use subscripts 2 and 1 to denote their y and x directions, respectively, and 5y is the Kronecker delta. Use 
the identity cos(fc/Xi) — cos(kfx' 1 ) — 2 cos (fc/ri/2) cos (kfX\). The ensemble and temporal averages of <&jj are 
2Fcos(fc/ri/2)cos(A;/Xi) [S i2 (Uj - C/j) + S j2 (V, - [//)], the trace of which is 4Fcos (fc/n/2) cos (k f Xi) (U 2 - U^). 
The X-space average of the first term in $y is 2F5i 2 cos (fc/ri/2) y J [J J (uj — u'-) dX 2 dX 3 ] cos (fc/Xi) dXi; inter- 
change i and j to obtain the second term in Whichever average is employed, this force introduces a term that 
has no small-scale spatial variation and is negligible in our scale analysis. 

Forced turbulence is temporally intermittent such that a space average, e.g., (da) M , does not obey dt (du} R = 0. 
The temporal intcrmittency observed by Borue and Orszag (1996) illustrates this fact; of particular relevance is the 
observation of repeated events characterized by accumulation of space-averaged energy in their mean flow (defined by a 
surface average in their calculation) , followed by a burst of transfer of energy from their mean flow to the space-averaged 

turbulent energy. Given the conditions mentioned below (p9j), one can time-average ( |7l| ) such that( dt (da) \ can 

be neglected; the time average has the effect of averaging the temporal intermittency. Now, apply to (|7l|) the r-sphere 
volume average ([43|). Consider the case in which the Reynolds number is large enough that r$ is in the inertial range, 
then we can neglect the term proportional to v in ( |7l| ) , and as shown in the preceding paragraph, any forcing term can 

be neglected. For a sufficiently long time average we have the approximation that § r (^(dm-^j ^ ~ — | ({^)s) T fs 

(recall that this is based on neglecting the time-derivative and viscous terms in (^) and the forcing because the 
forcing has no small-scale spatial variation and is therefore negligible in our scale analysis). A similar generalization 
of Kolmogorov's 4/5 law, namely, § r Da n ~ — | (e) E r, was obtained in Lindborg (1996), Frisch (1995), and Hill 
(1997) for the inertial range of homogeneous, anisotropic turbulence. 



Appendix B: Data 
B.l Inner Scales 



Inner scale was first defined by Obukhov (1949) as the r at which the asymptotic formulas for the inertial and 
viscous ranges are equal. Inner scales are more applicable in our scaling analysis by a factor of about 10 compared 

with rj. Inner scales for Dn and Dp/3, denoted t\\ and can be related to rj using (e) E = 15^^(9iUi) 2 ^ = 

(15/2) v (^(diup) 2 ^ , which is valid on the basis of local isotropy and incompressibility. Then, in — (2/3) 3 / 4 -^ = 
13r). Inner scales for Dm and D\pp, denoted £m and can be related to -q on the additional empirical basis that 

the derivative skewness, (^(diui) 3 ^ / ^(<9i?ii) 2 ^ , varies little from —0.5 over observed values of Reynolds number 

(Sreenivasan and Antonia, 1997; Belin et ai, 1997). Then = — 9. 677. The average of these four inner 

scales is (In + Ipp + 4u + Pipp) /4 = IO77. Even though the turbulence being studied need not be locally isotropic, 
we define I = IO77. This is the background of definition ([57[). 



B.2 Typical Data 



Data are needed for the investigations in Sec. 7.4. Let r/r, 1, and e be orthogonal unit vectors. Let subscript 
1 denote projection in the direction r/r, e.g., Dn = (r</r) (rj/r)Dij. Let subscript /3 denote projection in either 
the 1 or e directions (we need not distinguish which direction), e.g., Dppi is either (r^/r) -Dyfc or e,ej (rfc/r) Dyfc, 
but not e-iij (r^/r) Dijk. If a distinction need not be made as to the direction of projection, then subscript a is used; 
thus, D aa is either Dn or Dpp. No summation is implied by repeated Greek indices. A unit vector a in an arbitrary 
direction is a linear combination of the unit vectors r/r, 1, and e. Thus, if projections of a quantity have the same 
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order of magnitude and sign in all three directions r/r, 
has that order of magnitude. 

For the inertial range we use the formulas 



i, and e, then the projection in an arbitrary direction a also 



D aa = (e) 2 J\ 2/3 K aa (X,r,t) 
D ni = -{e) E rK in (X,r,t) , 



D 



1/3/3 



(z) E rK X pp ( X > r >*)' 



(Bl) 
(B2) 



The dimensionless coefficient functions, -Kqc, (X, r, t), ifm (X, r, i), and K\pp (X, r, i), are included to emphasize 
that our inertial-range data, like real data, need not be precisely homogeneous, locally isotropic, or stationary. The 
coefficient functions are assumed to be of the order of unity, and when differentiating the structure functions with 
respect to r j, the deri vati ves of the coefficient functions are assumed to be negligible compared to the derivative of 
r 2 / 3 in (Bl) and r in (B2). As motivation for this assumption, consider that for the case of local isotropy the above 
coefficient functions are constants between 2.7 and 0.26. The choice to scale with I = 10?y, rather than with rj, causes 
the coefficient functions to be of the order of unity. 



The slight effect of intermittency on the exponent 2/3 in (Bl) is not of concern here. Of more significance is 
the finding by Mydlarski and Warhaft (1996) of power-law ranges that are precursors to the inertial range. Their 
precursor power-law exponents are smaller than the 2/3 exponent of the inertial range, and the precursor exponents 
approach 2/3 as Reynolds number increases. Our scale analysis can be extended to apply for those weaker power 
laws; the accuracy of the scaling condition would be correspondingly weakened. 

Using I defined in (M), the viscous-range formulas for the scale analysis are 



D, 



act — x-~ I E 
6 



(efff/ 3 (r/£) 2 k aa (X,r,i), 



Dm = -(e) E l {r/tf fc m (X, r, t) , D X pp = -(e) E £ (r/lf k X pp (X, r, t) . 



(B3) 
(B4) 



where the dimensionless coefficient functions, k aa (X, r, t), km (X, r, i), and k\pp (X, r, t) are assumed to be of the 
order of unity, and when differentiating the structure functions with respect to n the derivatives of the coefficient 
functions are assumed to be negligible. In support of this assumption, note that for the case of local isotropy these 
coefficient functions are constants between 2.9 and 0.57. The choice to scale with I = 10?y, rather than with rj, causes 
the coefficient functions to be of the order of unity. 



D 



ill) 



and Dips have monotonic 



For r between the inertial and viscous ranges, the structure functions D a 
transitions between the asymptotic formulas (Bl), flB2|), and dB3|), (B4). Therefore, if a quantity is negligible on the 
basis of both (Bl), (B2), and (B3), (B4), then it is negligible for all r from within the inertial range to within the 
viscous range. 

Consider the projection of d r „Dij n in directions parallel to r, i.e., (fi/r) (rj/r)d rn Dij n , and perpendicular 
to r, e.g., eiejd, 

gives e^e. 



(d rn D, 
(n/r) (rj/r) (d rn D ijn ) 



Dijn 

For example, 

yn) = 9 r Dppi 



Note that neither projection commutes with the derivatives d ri , e.g., ei&j (d rn Dij n ) ^ 
if Dij n is locally isotropic, then differentiating the isotropic-tensor formula 
4 Dflfn, whereas <9 r „ (eiejDij n ) = d rn Dpp n + \Dpp\\ another example is: 



d r D 1 



li 



'/3/31, 

iDni - 



l -D 



0/31, 



i]n) — vr^j^ppn -f --Up/31, 

whereas <9 r „ [(^/r) (rj/r)D ijn ] = d rn D lln 

A i \ 

are 



d r D 1 



li 



-D 



in- 



I (e) E for the locally isotropic case. Since the 



By use of (B2), both e^ey (d rrl Dij n ) and (fi/r) (rj/r) d rn Dij n 

projections in all three directions are the same, the projection in the arbitrary direction a, i.e., aifiifl, 
equals — | 



for the locally isotropic case. 



J^^n^iju, also 

Although our data are not locally isotropic, assume that for our data 



di q>j df 



Dij n is of the order of — (e) E in the inertial range. For the viscous range and for the locally isotropic case, 



'1311 

Dijn — d r Dppi 



Dppi = fDppi = ^D ln and (n/r) (rj/r) d Tn D ljn = d r Dm + ~£>m - 7 -D/3/31 

(r/£) 2 for the viscous range. 



these formulas combined with (|B4|) imply that aia,jd rrl Dij n is of the order of — (e) E 



Assume that this is true of our data as well. 

The same manipulations apply to the projection of d Trl d r 



{d r + I) d r D u -r 
73/3 -;§r (-D/3/3 - Du). 



4 

73 



Dij . Differentiation of the isotropic-tensor formulas 



Dpp — Du) [in contrast, d rn d rn Du = (d r + -) d r Dn], and 
we find for the projection in an arbitrary 



By use of (Bl) and 



gives (n/r) (rj/r)d rn d rn Di 
ei&jdr n d rn D i: j = (d r + |) d r D 

direction a, that aiCijdr n dr n Dij is about 2 (e) 2 J 3 r -4 / 3 in an inertial range and is about 10 (e) E 
the viscous range. Assume that this is true of our data as well 



2/3 / _4/ 3 



. /2v in 



In effect, our definition of the inertial range includes (Bl), (B2), and that the projections in an arbitrary 
direction a of the two terms in ( f2~l| ) behave as 



cuajdr D 



13 n 

2uaia,jdr n dr n Dij ~ 0.2 (e) E (r/£)~ 4 ^ 3 . 



- e 



E ■ 



(B5) 
(B6) 
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where ~ means "is of the order of." Our definition of the viscous range includes (B3) and (B4), and that the 
projections in an arbitrary direction behave as 



a,iajd rn Dij n ~ - (e) E (r/£) , 
2vaiajd rn d rn Dij ~ (e) E . 



(B7) 
(B8) 



For both ranges we include the additional assumption that the off-diagonal components D a p (for a ^ 0) are not 
greater in magnitude than D aa . 

We assume that diagonal components of E^ (X, r, t) are of order (s) E ; that is, when projected on a arbitrary 
direction a, 



a,i&jEij ~ (e) E . 



(B9) 



In support of this assumption, recall that En = 4 (e) E on the basis of homogeneity, and, in the case of local isotropy, 
E aa =4{e) E /3. 

Data are needed for the diagonal elements aidfTij. Because ((p — p') (uj — u 'j)) E vanishes on the basis of 
local isotropy and because dx in operating on any average vanishes on the basis of homogeneity, it is assumed that 



ai&j (2 (p - p') (s i:j - s' i:j )) > ai&jdxi {{p - p') (uj - u'j)) 



(BIO) 



Of course, this is not true for the sum of diagonal components because of (|12|). It is likely that the right side of (B1C) 
is much smaller than the left side, but a more restrictive condition than ( p310|) is not needed. On the basis of DNS, 
Borue and Orszag (1996) show the cross spectrum of velocity and pressure gradient, where both velocity and pressure 
gradient are projected in their y-direction. For the inertial range, their data show that the corresponding diagonal 
component of T^- is proportional to (e) E r/L. 



More details for other flows would be welcome because TJy vanishes 



for locally isotropic turbulence (Hill, 1997). Therefore, its anisotropic behavior is of interest, 
result by Borue and Orszag (1996), it is assumed that our data in the inertial range obeys 



However, based on the 



(r/L). 



(Bll) 



The data by Alvelius and Johansson (2000) are consistent with (Bll). Using data from nearly homogeneous tur- 
bulence, Lindborg (1996) found that the single-point pressure strain correlation has a longitudinal component that 



(e) E /3 and a transverse component that is approximately 2(e) E /3. For homogeneous turbu- 
oo Lindborg's result corresponds to Tn — > —16 (e) E /3 and Tpp — > —8 (e) E /3; this agrees in 



is approximately —4 
lence in the limit r - 

order of magnitude with (Bll) evaluated at r = L. The first nonvanishing term of the Taylor series expansion of 



aiCij (2 (p — p') (sij — s'ij)) E is r 2 times the average of the product of pressure gradient and strain-rate gradient. This 
suggests that for the viscous range, 



.{r/lf{i/L). 



(B12) 



I. In the absence of further information, (B12) is assumed to be 



The form of (B12) is chosen to equal (Bll) at r 
valid. 

Finally, data are needed for Ti n and Tij n , which are defined in ( p3| ) and appear in the second term of ( |52[ ) . As 
described in Sec. 7.4, the second term in (|5|) is important for the case of large-scale structures in the flow. Assume 
that the experimenter chooses a flow that mitigates against large-scale structure; grid-generated turbulence is an 
example. In this case, it is assumed that 



< v 2 , and T ijn < v 3 , 



(B13) 



where v is defined in (p6[) . 

Relationships (JbTJ)- ( B 13 ) serve as exemplary data in the scale analysis. The fact that we have data only 
for projections in an arbitrary direction a means that we can investigate only the diagonal components of j2l]). 
The off-diagonal compon ents , whi ch are obtained by projection in two orthogonal directions, cannot be studied 
here. Relationships like (B1)-(B9) are most often associated with the assumption of local isotropy. However, like 
anemometry data, these relationships can be fulfilled for coefficient functions [as defined in (B1)-(JB4|)] of the order of 



unity without the specific restrictions of local isotropy being precisely fulfilled. For instance, for Dy, the restrictions 
for local isotropy are that its off-diagonal elements are zero, and that liijDij = eiejDij, and that Dn is related to 
Dpp by an incompressibility condition. In the scale analysis, such restrictions are not used; therefore, local isotropy 
is not assumed. 
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Appendix C: Homogeneity Implemented Using the Calculus of Local Homogeneity 

Although homogeneity is mentioned only briefly in this study it is useful to introduce it and to show how 
the calculus of local homogeneity produces the predictions of homogeneity for the case of homogeneous turbu- 
lence. Homogeneity is the approximation that ensemble averages do not depend on the position at which the 
average is obtained (Monin and Yaglom, 1975). That position being X, we implement this approximation by 
neglecting the result of dx n operating on any average. For example, in ( |2~l| ) \dx n dx n Dij = 0, dx n Fijn = 0, 
and for the average of (||) dx t ((p ~ p') (uj — u'j)) E — such that = — 2 ((p — p') (s^ — s^-)) „, from which 
we obtain Ta — because su — by incompressibility. The Taylor series of p(x, t) around point X is 
p(x, t) = p(X, t)+(x n — X n )dx n p(X,t) + ■■■. Upon averaging, homogeneity requires that dx n (p(X, t)) E = 0, etc., 
such that (p(x,t)) E = (j)(K,t)) E , and similarly (p(x',t)) E = (p(X, t)) E ; similarly, (e) E = (s') E . Within the average 
of ( ^o|) we have dx n dx n (p + p')e = 29x„^X n = 0, etc., such that En = 4(e) B . From homogeneity 

gives the incompressibility condition, d rn Dj n = 0. 

In Sec. 3.3 the spatial average is a volume average in X-space such that the equations [e.g., (|3^) and ([53])] do 
not contain dx n operating on an average. For those spatially averaged equations, homogeneity can be implemented 
by neglecting any average over the surface bounding the averaging volume of the surface-normal component of any 
vector. The basis for this implementation is that there are no net average fluxes in homogeneous turbulence. 



